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ABSTRACT 

Zonally-averaged autocovariances of the five-day height tendency for 
lags 0, 5, , . , , 25 days are computed at 25 rib around latitudes 55N, 65N, 
and 75 R during the period 5 October 1958 through 29 March 1959. Using 
the convolution theorem for each space-wave number, it is possible to 
determine the zonally-averaged spectral density of the height-tendency 
field at each of these latitudes in terms of tine waves having periods 
from 10 to 50 days. The results indicate that all but about 1 0,1 of the 
variance of the height-change field is accounted for by space waves one 
and two. For space wave number two, time waves number 3 end 4 (per 50 
days) give a significant contribution to the spectrum at the 90/5 confi- 
dence level. However, for space wave number one, time wave number one 
gives the most significant contribution to the spectrum. It is inferred 
from this, that cases of "explosive" warning, which may occur as early 
as January are correlated with intensification and/or movement associ- 
ated with space wave number one, at periods of approximately 100 days. 

The author wishes to express his sincere appreciation to Professor 
Frank L, Martin, Department of Meteorology, U. S, Naval Postgraduate 
School, for his suggestion of the topic and his encouragement in the 
preparation of this study. 
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1 . Introduction 



During the past ten years, meteorologists have become aware of the 
winter-time phenomenon known as "explosive" warming, which occurs in the 
Arctic and sub-Arctic stratosphere. During this same decade, the increased 
availability of higher-level soundings has made accessible to the meteor- 
ologist sufficient information to analyze the structure and causes of 
these unusual occurrences. 

The work of the McC-ill University group, primarily under P. K. Hare 
[i960] and B. W. Boville [i960] has made it clear that there is a succes- 
sion of cold troughs and warm ridges at the 25 -mb level which eventually 
culminates in the development of an intense warm anticyclone in polar 
latitudes, northern Canada is a preferred location for the occurrence of 
this phenomenon. This intense anticyclogenesis may be explained JjBoville, 
I960] in terms of baroclinic instability in the polar night westerlies. 
Actually the development of the intense anticyclone represents a stage in 
the weakening of the polar night westerlies [Julian, 1959] which are re- 
placed in spring and summer by the easterly wind regime that finally en- 
velops the entire summer stratosphere. Boville and his associates [I96l] 
have prepared an atlas of the 25-mb circulation for the six-month period 
5 October 1958 through 29 March 1959. In these studies, Boville et al. 
list the Fourier analyses of the 25-mb contour field at five-day intervals 
for latitudes 55N, 65N, and 75N. For each of these latitude circles, they 



describe the Fourier analyses in the form 

9 



z = z +Jc n co5(nA-(j) n ) 

n=1 



( 1 ) 



where the values of C^, and are listed in the atlas through space-wave 
number nine. Here "2, is the zonal mean contour height around the latitude 
circle while and (J)^ are the amplitude and phase angle, respectively, 
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of the nth harmonic (n=1 , ,9). 

Spectral analysis has become a useful diagnostic tool of meteorolo- 
gists, primarily for the purpose of detecting preferred periodicities in 
the eddies or perturbations under consideration. The number and variety 
of such studies in certain areas of statistical meteorology are too great 
to attempt to describe them here. Suffice it to say that the basic methods 
were largely developed by J. U. Tukey, and a concise summary of their ap- 
plications to Meteorology can be found in the excellent textbook by 
Panofsky and Brier £l958]. 

Actually, Boville [i960] has published the spectral analysis of the 
2 5 -mb contour height field at latitude 60N for the month of October 1958. 
In that analysis, one of the conclusions reached was that most of the 
variance in the contour height field at 25 mb is explained by space-wave 
numbers one and two. A similar conclusion is reached in this paper in 
regard to the variance of the height tendency. 

The objective of the present paper may be understood from the local- 
time differentiation of equation (l ) 




n=l 



where i”*- <3z :/dt- dz/dt . From equation (2), it may be 
seen that time-changes in the perturbation field may be ascribed to terms 
involving C^, the local-time change in the harmonic wave amplitude, and 
(J) n , the local-time change of the phase angle. These two types of terms 
in equation (2) may be described, as the int ens if i cat i on and movement con- 
tributions, respectively, for the nth wave. A time-spectral analysis of 
d^/dt, zonally averaged with respect to longitude, was then performed 
using the five-day rates of change and (J) n . (For the mathematical de- 
velopment see section 2. ) An advantage of the spectral analysis of the 
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height-tendency field over that of the height field alone is that the 



former can indicate periodicities due to intensification and to movement 
of the harmonic waves, whereas these connections are not so clear in the 
latter case. 



3 



2. Mathematical development 



The autocovariance of Z* at lag X is given [Blackman and Tukey, 

1958] by the formula 

/r/2 

R(X)= 6 i / ilt) 2 .*(t+T)dt. (3) 

7^->0O 7- J 

J-T/Z 

where Z* is the zonal height perturbation, and includes the contribution 
from all space waves around the latitude circle, T is the averaging per- 
iod over which the autocovariance is determined, X is the lag, and t is 
an interval of time, also measured in five-day units. Using the truncated 
time series available for this study, the time integration may be replaced 
by a summation over < Tl-K equally-spaced time intervals, at lag X = KA t . 
Then the integral (3) may be rewritten in the form 

JM ; . o^K-M, (4) 

u t=i 

where R (K ) is the zonally-averaged time autocovariance. In equation (4), 

A represents the longitude at any point within a latitude circle. Here, K 
is me number of lags in terms of intervals of At = 5 days , Yl is the num- 
ber of five-day height changes available in the time series, and M is the 
maximum number of time lags used in this study. Obtaining the five-day 
perturbation height tendency + at time t + X by an equation 

analogous to (2), we may write, 



2T TtK 



R(k)=_l( ±)l 






3 

) =JJi cos (n A- tfv) + a, ift sm(n A - %)l (5 ) 

n~i 



where Q>, and (j//, are analogous to Cn and (f> n of equation (l). Hence, 
from equations (2), (4), and (5), the zonally-averaged time autocovariance 
for space wave n has the form 
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( 6 ) 



fflds. 



Rr.iKj= x \ f[ . 4 G4 fa hi co: f ^-4) 
%~k / ; ! 7 j 






In equation (6) the subscript n indicates that R^K) is independent of all 
space waves other than the nth. In the derivation of (6), considerable 
simplification has been made possible as a result of orthogonality prop- 
erties, that is by deleting terms which vanish upon integration around a 
complete latitude circle. A proof of equation (6) is given Appendix I. 
Since discrete data and a finite time series have been used,' the maximum 
lag Mat limits the maximum resolvable period to be Xa« = 2M At , 
which may be called the fundamental period . In this study the maximum 
number of lags has been taken as M = 5 so that the fundamental period will 
be 50 days. If K is allowed to range over the values 0, 1, 2, . .., M, 
various values of R^K ) can be generated for the different space-wave 
numbers. Blackman and Tukey [l95Sj recommend, when using a truncated 
time series, that the maximum number of lags M should not exceed approx- 
imately 1 0 % of the length of the record. Since the total available number 
of five-day mean values of c) Z. /() ~t was 35> the maximum number of lags 
was taken to be five. This lag number was considered to be in reasonable 
accord with the recommendation of Blackman and Tukey, and yet sufficiently 
large to give resolution for the long-period waves which occur in this 
meteorological regime. 

The definition of the spectral density at time -wave number N [Kahn, 
1957^ is the average contribution to the variance from a unit wave number 
interval centered at wave number N. Since the maximum period under con- 
sideration is 50 days, the concept of a time-wave number N is to be inter- 
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preted in terms of frequencies as follows: 



1 wave per 50 days for N = 1 , 



N waves per 50 days for wave number N, N ^ 5. 
Thus the frequency may be expressed as 



= _£L_ , P'1, (v) 

zmt 7 

and spans the frequency-interval (N - -^)/2MAt to (N + fTAt. The 
corresponding period for this time-wave is the reciprocal of . The 
highest frequency resolvable when using such discrete data at finite 
intervals occurs for N = M and is 





is often called the folding or Nyquist frequency, and in the present 
study is one cycle (l 0 days ) . 

Panofsky and Brier [l958j give the following equation for raw spec- 
tral density estimates for time-wave number N: 



5 N = 2 ji [R(o) + (-1)“ R(m)] 




— M-l.(s) 



In practice, S N is generally considered to be distributed uniformly over 
the unit wave-number interval N - •§• to N + -i,-, except at N = 0 and N = M. 

For the special cases N = 0 and N = M, the corresponding estimates 
S Q and S M follow from equation (8) multipled by one -he. If because only 
one-half of a wave-number interval is available at N = 0 and N = M. 

These unsmoothed raw spectral estimates may be smoothed by "harming" 
[Blackman and Tukey, 195&], in the following manner 




r 



S' + ^ s + 

mm Jm < 






s 

U/JtI 



(9) 
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except at the extreme wave numbers N = 0 and N = M where the average is 
modified to 




( 10 ) 



These smoothed estimates may be considered to be centered at wave number 



Note that the integrated value of all the spectral estimates for a 
particular space-wave n gives the total time variance, for this wave* In 
other words, by definition, the relationship 



holds, although carrying the integration only to N = 5 neglects some of 
the higher frequency spectral contributions. In addition, "aliasing" or 
folding of high frequency energy back upon the Nyquist frequency has been 
neglected. However the spectral densities at higher space-wave numbers 
(table 2) suggest that this is not important in the present problem. 

This has been inferred since the contributions to the time variance for 
these wave numbers appears to act as "noise", with no preference for the 
higher time-wave numbers. 



N, or frequency 




equation (7). 




( 11 ) 
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3 . Computational procedures and interpretations 

As mentioned in section 2, the maximum number of time lags chosen 
was M = 5. This choice gives a resolution of the spectrum of bi/bt 
to a maximum frequency of l/lO cycle (day) ' . 

All computational processes were performed on the CDC-1604 computer. 
All programs devised for this study we re written by the author or modified 
from those available in the U. S. Naval Postgraduate School computer cen- 
ter library. The three latitudes chosen for analysis were 55N, 65N, and 
75N. Table 1 , a, b, c gives the zonally-averaged autocovariances for each 
space-wave number computed by equation (6) at each of the three latitude 
circles. The smoothed density spectra for each space-wave number were 
then computed according to equations (8), ( 9 ), and (lO) and the results 
displayed in table 2, a, b, c. 

Note that in table 2, a, b, c periodicities are presented down the 
left-hand column, while the space-wave numbers n = 1, 9 appear along 

the top of the table. From the rows denoted "P. R. Var.", it may be seen 
that the space waves n = 1 and n = 2 together account for approximately 
90$ of the variance of the zonally-averaged five-day height changes. 

While the percent reduction of the variance of dz./d"t accounted for by 
space wave 2 exceeds that of wave 1, the contribution of the latter in- 
creases as one proceeds from latitude 55 N to 75N. The contribution by 
waves 3 and 4 accounts for less than 10% at latitudes 55 N and 65N, and 
for less than 5% at latitude 75 N. Thus, in this study only the time- 
spectra of space waves one and two have been analyzed for significant 
periodicities . 

For space waves n = 1 and n = 2, fiducial confidence limits at the 
10% and 90% significance level have been taken from Blackman and Tukey 
[1953]. These limits are based upon the t/i distribution, where 
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(a) 55N 



MG SPACE HAVE NUMBERS 

DAYS 





1 


2 


3 


4 


5 


6 


7 


8 . 


9 


0 


42.18 


218.66 


11 .80 


13.43 


5.35 


1.64 


2.01 


2.04 


.74 


5 


8.86 


-31 .97 


2.63 


-1 .28 


.31 


.14 


.01 


.14 


.04 


10 


2.43 


-31.53 


4.22 


.75 


.26 


-.02 


.16 


.19 


.04 


15 


-3.37 


6.46 


1.71 


2.55 


.79 


-.06 


.32 


.45 


.04 


20 


.86 


5.03 


3.25 


.44 


-.72 


.14 


.40 


.37 


-.03 


25 


~8 .46 


-11 .47 


2.22 


1.29 


.36 


.30 


.32 


.54 


.14 
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(b) 65N 












0 


80.78 


336.38 


15.31 


27.36 


3.11 


.97 


1.37 


1 .14 


.27 


5 


23.62 


-139.61 


1.25 


5.23 


.53 


.12 


-.01 


.07 


.01 


10 


-1.55 


25.27 


-2.33 


8.17 


.07 


-.02 


.58 


.09 


-.02 


15 


-4.05 


14.77 


.50 


5.78 


.04 


-.01 


.02 


.18 


-.02 


20 


1.33 


-13.28 


.58 


-.28 


.12 


-.03 


,23 


.24 


.07 


25 


-17.53 


-1 1 .36 


1 .18 


3.94 


-.10 


.00 


.08 


.11 


.03 










(c) 75N 












0 


103.59 


139.47 


6.21 


2.69 


1 .39 


1 .07 


.34 


.36 


.23 


5 


28.08 


-31 .36 


.97 


.44 


-.04 


.22 


.02 


.10 


.04 


10 


7.89 


4.73 


-1.02 


.17 


-.01 


-.03 


-.00 


.03 


.06 


15 


-11.79 


-1 .13 


-.10 


.79 


.06 


.09 


.03 


.02 


.02 


20 


-4.05 


42.43 


.19 


-.23 


.25 


.05 


.11 


.04 


-.03 


25 


-11 .88 


-46.88 


.41 


-.26 


.01 


.20 


.09 


.04 


-.00 



Table 1. Zonally averaged autocovariances of bi/ bt for lags of 0, 5, 

. .., 25 days at latitudes (a) 55^, (b) 65N; (c) 75W. The units 
are in terms of 10 4 (feet/5 days) 2 . 
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(a) 55N 



FREQUENCY 






SPACE 


WAVE 


NUMBERS 








CYCLES/DAY 


1 


2 


3 


4 


5 


6 


7 


8 


9 


0/50 


9.28 


19.34 


2.89 


1 .81 


.98 


.25 


.30 


.34 


.13 


1/50 


10.06 


23.81 


2.36 


1 .76 


1 .00 


.31 


.29 


.31 


.14 


2/50 


9.30 


46.17 


1 .61 


2.11 


.95 


.36 


.34 


.34 


.14 


3/50 


6.57 


54.77 


1 .40 


2.96 


1.06 


.31 


.41 


.41 


.14 


4/50 


4.71 


44.96 


1.51 


2.92 


1 .01 


.24 


.36 


.36 


.13 


5/50 


4.05 


34. .94. 


1 .61 


2.41 


.79 


.21 


.30 


.29 


.12 


P. R. Var. 


14.2 


77.0 


3.5 


4.5 


1.9 


0.7 


0.6 


0.7 


0.3 










(b) 65N 










0/50 


18.69 


19.95 


2.17 


6.91 


.62 


.18 


.28 


.19 


.03 


1/50 


20.53 


25.75 


2.79 


5.78 


.67 


.20 


.23 


.18 


.04 


2/50 


19.59 


44.51 


3.65 


3.78 


.66 


.21 


.15 


.20 


.06 


3/50 


13.41 


79.03 


3.48 


3.90 


.55 


.18 


.16 


.22 


.06 


4/50 


7.44 


95.69 


2.33 


4.52 


.41 


.14 


.25 


.19 


.04 


5/50 _ _ 


4.88 


90.79 


1 .60 


4.13 


. *33 


.11 


.29 


.15 


.03 


P. R. Var. 


16.0 


75.0 


3.1 


5.2 


0.6 


0.2 


0.2 


0.2 


0.1 










(c) 75N 










0/50 


24.56 


13.15 


.98 


.60 


.18 


.21 


.05 


.08 


.06 


1/50 


26.90 


14.79 


1 .27 


.53 


.21 


.22 


.05 


.08 


.05 


2/50 


23.46 


24.01 


1 .58 


.46 


.27 


.24 


.07 


.08 


.03 


3/50 


14.54 


30.77 


1.34 


.54 


.29 


.21 


.07 


.06 


.03 


4/50 


10.50 


27.46 


.81 


.46 


.24 


.14 


.05 


.04 


.04 


5/5Q 


_9_.57 . 


24.97 


.53 


.29 


.20 


.10 


.04 


.03 


.03 


P. R. Var. 


42.0 


53.0 


2.6 


1 .1 


0.5 


.4 


0.1 


0.1 


0.1 



Table 2. Zonally-averaged spectral densities for space-wave number n=1 
...,9, and wave frequency N=0,l/50, . . .,5/50 cycles/day. Sec 
tions (a), (b), and (c) are for latitudes 55N> 65N, and 75N. 
The units are in terms of l0 4 (feet/5 davs) 2 . 
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/= (2 9(- M/Z)/n is the nunber of degrees 



is the nunber of degrees of freedom. In this in- 



vestigation we have *77 = 35 and M = 5, so that f ^ 



14. For these values. 



the 10,9 and 90'" fiducial limits interpolated from those given by Blackman 



Figure 1 illustrates the time-spectral analysis at 65N, for space 
waves n = 1 and n = 2 (the latter shown by the solid line). From this 
diagram it may be seen that for n = 2, the nean variance averaged over all 



structed relative to this mean. Thus it is evident that for n = 2, waves 
of period 10 days and 12.5 days give significantly large contributions at 
a 90 % confidence level. On the other hand, at 65N space wave 2 at a per- 
iod of 50 days, gives significantly small contributions to the variance. 

The same type of analysis for n = 1 at latitude 65 W shows that the 
mean variance for space wave n = 1, averaged over all -periods is 13.5, in 
units of l0 4 (ft/5 days) , so that a time-wave of period 50 days gives a 
significant contribution at a 90%> confidence level. It should be recalled 
that a period of 50 days in figure 1 actually spans a range of periods 
from 33.3 to 100 days, and such waves could include those associated with 
explosive wanning. This is so, since strictly speaking, the largest value 
of bC | /c) t which contributes to the variance of bz/^t could occur at 
a time which is either three-fourths or one and three-fourths of a period 
after the onset of the autumnal cooling season in polar latitudes. 

In table 3 the mean spectral densities at latitudes 55N and 75N, 
averaged over all frequencies in table 2, are given for space waves n = 1 
and n = 2. 




periods is 60.1 x l0**(ft/5 days )^. The fiducial limits have been con- 
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0 1/50 2/50 3/50 4/50 1/10 FREQUENCY (CYCLES/DAY) 

50 25 lb.l 12.5 10 PERIOD (DAYS) 

Figure 1. Smoothed time spectral densities in units of l0^(feet/5 days) , at 65W» The solid line is for 
space wave n = 2, while the dashed line is for space wave n = 1 • 



LAT 


SPACE ■ 


- WAVE 




1 


2 


55 


7.4 


40.4 


75 


18.5 


23.2 



4 2 

Table 3. Mean spectral densities in units of 10 (feet/5 days) at lati- 
tudes 55N and 75N for space-waves n = 1 and n = 2. 

Using the fiducial confidence limits for the 2/^ distribution pre- 
viously mentioned (0.55 and 1.51, respectively), it is apparent from table 
2 (a, c} that there exists no significant preferred peak in the spectral 
densities at the 90/ limit for either latitude 55N or 75N, 

For space wave n = 1, at each of the latitude circles of table 2, 
there is a steady decrease of spectral density, S N , with increasing time- 
wave H ^ 1 . 

For the case of isotropic turbulence, MacCready [l 953] has tested the 
applicability of a spectral function of form 



5 N cC N 



where N is the tine-wave number. This test is essentially based upon a 
theory of isotropic turbulence set forth by Heisenberg |l 948j in which 
p = 5/3. For n = 1, and using the results of table 2, computation of the 
values of the exponent p gave the following: 



55N, 


p = 0.57 


65N, 


p = 0.91 


75N, 


p = 0.64 



Chiu [l96oJ has obtained somewhat similar results for the exponent p using 
the spectrum of the zonal and meridional components of the wind at Belmar 
and Cocoa, and finds values of p of the order of 1 . 

On the other hand, the results of this study for space wave n = 2 do 
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not show the decrease of spectral density \;ith increasing time-wave num- 
ber. Hence for n = 2, the peak spectral density which occurs in each case 
of table 2 at either time-wave 3 or 4 is even more likely to be signifi- 
cant of the average state of synoptic events. 

Space-wave n = 1 , with its maximum spectral contribution at time -wave 
N = 1 is capable of explaining those perturbations of longer period asso- 
ciated with "explosive" warming. 
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4. Conclusions 

(l ) Space waves 1 and 2 contribute 90£> of the variance of the height- 
tendency field at latitudes 55*’, 65N, and 75N. 

(2) The chief contributory to the spectrum of hz/c)t occur for 
tine -wav ey N = 3 or N = 4 (periods 12.5 and 16.7 days), which have the same 
general periodicity as the average perturbation period at a fixed station 
jjsee, for example, Austin and Krawitz, 195&J. 

(3 ) The chief contributor to the time variance for space wave n = 1 
occurs for time wave II = 1, or possibly N = -y, which corresponds to periods 
of the order of 50 days or longer. Such periods may be sufficiently long 
to account for the "explosive" wanning which sometimes occurs as early as 
January . 

(4) For space wave n = 1, the time spectral density was distributed 
with time wave number according to 

r n 1 - a7 ' 

which compares reasonably with the results obtained in other spectral 
studies Jc.f. Chiu, 196 oJ. 
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APPENDIX I 



DERIVATION 0? THE ZONALLY-AVERAGED TIME AUTOCOVARIANCE 



The equation for the zonally-averaged tine autocovariance of the 
height-tendency field at lag t nay be expressed in the form shown in 
equation (4 ) „ 

/ZW ttn 



R(w= l ( j _ ) fz*(t) z*(t+x)l6A, 

nfh-h/J- J 



( 4 ) 



J o t=1 

where T 1 is the maximum number of discrete values in the sample, and R(k) 
is the zonal ly-avcraged time autocovariance. 

When equation ( 4 ) is expanded using equations (2) and (5), and upon 
interchanging the order of summation, equation (12) is obtained 




Cn j L [sin (2 nA-i)n-y + S,N ( llVr &„)] 



+ 



(tn%C n fiilN f 2 nA - (f) n - l)J ) - S'f, I 



i! + 



£ 



(4 L, UJ, (cu.yll, - (bo) - cos (2 n A — c|) n - V)J ( . )] A/' 



(12) 



where the orthogonality property has also been used. Equation (12) repre- 
sents a double integral; first, with respect to longitude, second, with 
respect to time. The integral with respect to longitude gives the zonal 
average. On integrating around the latitude circle from 0 to 27T , all 
cosine and sine terms involving the argument ( 2 . nA — (j) h — (jj,^ 
integrate to zero. Therefore, the result obtained is 
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o, 









[oJainlO.Mi] &H'(l|'n-$n) | 



(6) 



Note that Pn(K) is the sonally-averaged autocovariance at a tiir.e-lag 
of KAt days, for the nth space-wave; no other space-wave appears in 
the expression for R n (K). The interval between successive data is 

r 

such that A 5 days . 
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